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The Very High Energy Radiation Telescope Array (VERITAS) is a system of four imaging Cherenkov tele- 
scopes currently under construction at Kitt Peak, Arizona, USA. The first telescope has been in operation at 
the Mt. Hopkins basecamp since January 2005. We present here detailed Monte Carlo simulations of the tele- 
scope response to extensive air showers. The energy threshold for this stand-alone telescope is calculated to be 
150 GeV at trigger level, the 7-ray trigger rate is 22 7's/min. Image parameter distributions, and the quahty of 
gamma-hadron discrimination are calculated and show good agreement with distributions from observations of 
background cosmic rays and high-energy gamma-rays from the Crab Nebula and Markarian 421. The energy 
spectrum of the Crab is reconstructed as (3.26 ±0.9) • 10"^ . £;-(2.6±o.3) m-^s-iTeV^i. 

1. Introduction 

The first of the four VERITAS telescopes has been in routine operation since January 2005 at the basecamp of 
the Fred Lawrence Whipple Observatory at Mt. Hopkins, Arizona. The site is 1275 m above sea level. Several 
TeV gamma ray sources have been detected 1 1 1 and the construction of the other telescopes is underway. A 
detailed description of the first VERITAS telescope, called Telescope 1 in the following, can be found in (T], a 
fuU description of the VERITAS project in |2|. 

2. Data 

For the following comparison of data with Monte Carlo simulations, gamma-ray candidates are extracted from 
observations of the Crab Nebula (3.9 h on source) and Mrk 421 (4.2h on source) in March and April 2005. 
The data is taken in ON-OFF mode at elevations above 60° and in good weather conditions. Images are 
parameterized using standard second moment analysis |3|. Image pixels are selected by a two level cut on 
the charge per pixel. Picture and boundary thresholds are pedestal plus four respective two times the pedestal 
variance. The cuts on image parameters to select gamma-rays are a < 10°, length/size < 0.002, length 
> 0.14°, 0.067° < width < 0.1°, and 0.5° < distance < 1°. Alpha is the orientation angle of the image ellipse 
relative to the direction of the source position. Length and width describe the shape of the image, and distance 
the location of the image centroid with respect to the source position. Only events with at least five image 
pixels are taken into account. These cuts results in a combined significance of about ll.la and roughly 900 
gamma-ray candidates for the observations of the Crab Nebula and Mrk 421. The Crab data set is statistically 
limited due to the fact that the observation started very late in the season. 

3. Monte Carlo simulations 



A complete chain of Monte Carlo simulations has been developed. It consists of air shower simulations with 
CORSIKA j4l and a detailed simulation of the telescope response 0. CORSIKA version 6.20 is used with 
the hadronic interaction models QGSJet for primary energies above 500 GeV and FLUKA below. Simulated 
hydrogen, helium, and oxygen nuclei induced air showers are produced in an energy range from 10 GeV (150 
GeV for oxygen) to 50 TeV at different elevations. Spectral indices are taken from fits to balloon measurements 
l6l . The shower cores are distributed randomly on a circular area with a radius of 1000m around the telescope. 
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The isotropic distribution of the cosmic ray incident angles is simulated by randomizing the shower directions 
in a cone of radius 3.5° around the pointing direction of the telescope. Primary gamma-rays are simulated 
with a Crab-like spectrum in an energy range from 10 GeV to 10 TeV on a circular area at the ground with 
radius 450m. Scatter areas are chosen according to the lateral Cherenkov light distributions of hadronic and 
gamma-ray showers. 

Measurements of the atmospheric properties at the current site of Telescope 1 and the final one at Kitt Peak are 
currently in progress. The calculations described here use the U.S. standard atmosphere, which does not always 
reflect the properties of the atmosphere in southern Arizona. Atmospheric extinction values are estimated by 
MODTRAN 4 1 7 1 assuming 50 km visibility at 550 nm at ground level. The photo-electron rate per PMT is 
measured to be 100-200 MHz. This corresponds to a night sky background rate of 2.8 x 10^^ photons • mT^ 
s^^ sr^^, which is used in the simulations. 

The telescope simulations consist of two parts, the propagation of Cherenkov photons through the optical 
system and the response of the camera and electronics. The geometrical properties of the optical system 
is fully implemented in the simulation. Misalignment of the mirrors and their surface roughness is taken into 
account. The camera configuration is that of April 2005, i.e. a 499 pixel camera with 1 1/8" phototubes without 
light cones. The field of view of the camera is 3.5°. The response of the PMTs to single photons has been 
measured; the single photo-electron pulse has a rise time of 3.3 ns and a width of 6.5 ns 1 1 1. In the simulations 
a signal in a PMT is created by adding up single photo-electron pulses with appropriate amplitude and time 
jitters applied. Electronic noise and all efficiencies, including mirror reflectivities, geometrical, quantum, and 
collection efficiencies, and losses due to signal transmission have been modeled. The pulses are digitized 
into 2ns samples with a trace length of 24 samples reflecting the properties of the FADC system. The trigger 
simulation utilizes a simplified model of the constant fraction discriminator, which is the first stage of the 
VERITAS multi-level trigger, and a full realization of the pattern trigger, requiring three adjacent pixels above 
threshold in a time window of 5 ns. The currently used trigger threshold of 70 mV corresponds to about 6.7 
photoelectrons. The output of the telescope simulations, i.e. FADC traces for all PMTs, are written to disk in 
the VERITAS raw data format. The analysis steps, which include pedestal calculation, image cleaning, image 
parametrisation, and source reconstruction, are exactly the same for simulated and real data. 

4. Comparison of Data with Monte Carlo simulation 

4.1 Trigger rate 

The raw trigger rate of Telescope 1 with trigger conditions as described above is between 140 and 150 Hz, 
mainly depending on elevation. Dead time losses due to readout are ^10%, the corrected trigger rate is conse- 
quentially ^150-160 Hz. The simulation of the cosmic ray background results in trigger rates of 101 Hz from 
air showers induced by protons, 26 Hz from helium nuclei, and ^5 Hz from nuclei of the CNO group for an 
elevation of 70°. The Monte Carlo calculations reproduce the observed trigger rate with an accuracy of ^15%, 
which is acceptable, taking into account an estimated uncertainty in the cosmic ray fluxes of ^^25% and various 
systematic uncertainties in the modeling of the telescope. 

4.2 Image parameter distribution 

Second moment image parameters are the standard tool in analyzing data from lACTs. The suppression of 
the hadronic and muonic background as well as the estimation of the primary energy depends strongly on the 
correct predictions of the image parameters from Monte Carlo calculations. Fig^shows a comparison of image 
parameter distributions from real data and simulations. The alpha, distance, length, and width distributions 
show very good agreement and proof the validity of the simulations. 

The described cuts on image parameters are chosen to extract a sample of gamma ray candidates with very little 
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Figure 1. Comparison of image parameter distributions from observation data (closed symbols) and IVIonte Carlo calcula- 
tions (open symbols). 



contamination from background events. About 10% of all gamma ray showers are accepted with these cuts, 
rejecting ~99.9% of all cosmic ray showers. These hard cuts are necessary for data from a single telescope to 
reject the large background due to local muons. (Telescope aiTays are insensitive to such single muons and can 
be run and analyzed with looser cuts.) 

4.3 Performance 

The effective area of Telescope 1 after cuts is shown in Fig|2l(left). Most of the high-energy events fail the 
hard 7-selection cuts due to the cut on distance. This is different for the application of extended supercuts fsl, 
which do not result in a similar good rejection of background events, but take the dependence of the image 
parameters on image size into account. The maximum effective area after extended supercuts is ~ 2.8 • 10^ 
m^. The measured gamma ray rate from the observation of the Crab Nebula is 2.1 ± 0.2 7's/min. Using the 
Crab spectra reported in |9| results in a simulated gamma ray rate from the Crab Nebula of 2.2 7's/min. The 
energy threshold of a Cherenkov telescope is conventionally defined as the position of the peak of the energy 
spectrum of the source convoluted with the effective area curve of the detector. According to this definition, the 



4 



Gemot Maier et ah 




log,„ energy [TeV] log^^, energy [TeV] 



Figure 2. Left: Effective area of Telescope 1 after 7-selection cuts (filled circles) and extended supercuts (filled triangles). 
Right: Energy spectrum of the Crab Nebula compared with earlier results from Whipple 1 9 1 and HEGRA 1 10). 



threshold is 150 GeV at trigger level, 160 GeV after applying extended super cuts, and 370 GeV after applying 
the hard cuts described earlier. The right-hand side of Fig|2] shows the reconstructed energy spectrum of the 
Crab Nebula after applying extended supercuts to the data set. A power law fit to the data points results in a 
spectral index of 2.6 ± 0.3 and a flux constant of (3.26 ± 0.9) • 10^^ m~^s~^TeV^^ (statistical errors only). 
This agrees well with earlier measurements by other telescopes. 

The agreement of the Monte Carlo simulations with observational data shows that we have a good understand- 
ing of the first VERITAS telescope and that the design performance is being met. The final configuration of 
VERITAS with four telescopes will have a significant higher sensitivity and an energy threshold of 1 10 GeV 

CD. 
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